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ABSTRACT 

We report the first measurement of genus curves for the two-dimensional mass 
map in the neighborhood of rich, X-ray luminous galaxy cluster MS1054-03 at 
z = 0.83, reconstructed from weak lensing data obtained by Suprime-Cam on the 
prime focus of 8.2m Subaru telescope 1 . We find that the genus curve measured in 
the whole survey field deviates from that expected from a random Gaussian field. 

We show that this non-Gaussianity is induced by the rich cluster in this region, 
and that the genus curve for the region without the cluster is consistent with the 
prediction for a random Gaussian field. The results demonstrate the feasibility 
of the genus statistics to examine the non-Gaussianity due to the large scale 
structures and to probe the statistical properties of the large scale clustering. 

Subject headings: cosmology: observations - gravitational lensing - large-scale 
structure of universe - galaxies: clusters: individual(MS1054-03) 


1. Introduction 

The origin and the evolution of the large scale structures attract much attention and the 
understanding about the statistical nature of the mass distributions in the universe is one of 
the important issues to be addressed in the observational cosmology. The inflationary theory 
gives us a definite prediction on the statistical properties, namely the random Gaussian 


1 Based on data collected at Subaru Telescope, which is operated by the National Astronomical Observa¬ 

tory of Japan. 
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nature of the initial density fluctuation. The density fluctuation evolves by gravitational 
instability to have the non-linear distribution, i.e. non-Gaussianity, at later stage whose 
details may depend on the cosmological parameters as well as the nature of dark matter. 
Thus the measurement of statistical properties of the density distribution is very important 
for the understanding of the formation and evolution of large scale structures. 

The topology of large scale structures was qualitatively described by Einasto et al. 
(1984) and its quantitative study was initiated by Gott, Melott & Dickinson (1986), Hamil¬ 
ton, Gott & Weinberg (1986), and Gott, Weinberg & Melott (1987), and has been studied 
by many authors on the basis of various observational data. For example, Coles & Bar- 
row (1987), Melott et al. (1989), and Gott et al. (1990) investigated the genus or the 
Euler-Poincare characteristic statistics of the topology in two-dimensional fields of galaxy 
distribution. Further, the method has been applied for various kinds of observational sam¬ 
ples, the angular distributions of galaxies on the sky (Coles & Plionis 1991; Gott, Mao, Park 
& Lahav 1992), the galaxy distributions in slices of the universe (Park et al. 1992; Colley 
1997) and the temperature fluctuation fields of the Cosmic Microwave Background (Coles 
1988; Gott et al. 1990; Smoot et al. 1994; Kogut et al. 1996; Colley, Gott & Park 1996; 
Park et al. 1998). Recently the topology of the galaxy distributions in 2 degree Field Galaxy 
Redshift Survey and Sloan Digital Sky Survey has been studied by Hoyle et al. (2002) and 
Hoyle et al. (2002). These studies show that the topology of the two-dimensional galaxy 
distribution is consistent with the random phase Gaussian distribution with a possible weak 
’meat-ball’ topology where the topology is called ’meat-ball’ when it looks like collection of 
isolated islands. More recently, Park, Gott & Choi (2001) calculated the two-dimensional 
genus statistics for the contour maps of the galaxy distribution in the Hubble Deep Field. 
Using the red-shift distribution of galaxies measured by the method of photometric red-shift, 
they studied the evolution effect in the topology. However, the genus curves was consistent 
with random Gaussian field, and they were not able to discover any meaningful red-shift 
dependence. 

Previous topological studies using galaxy distributions suffer from the biasing effect 
between ’visible’ galaxies and underlying mass distributions and it is not clear whether 
the galaxy distribution traces mass distribution in the universe precisely. Thus it would 
be very useful to develop a method to measure mass distributions in the universe directly. 
We shall make use of weak gravitational lensing techniques for this purpose. Weak lensing 
directly measures the projected mass distributions from the observed image distortions of 
background galaxies and thus can be used to compare the observational results directly with 
the theoretical and/or numerical predictions. The reliability of the weak lensing techniques 
for measuring the mass distributions on large scales has been widely recognized, for example, 
by cosmic shear measurements (van Waerbeke et al. 2000; Wittman et al. 2000; Bacon, 
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Refregier & Ellis 2000; Kaiser, Wilson & Luppino 2000; Maoli et al. 2001; Rhodes, 
Regregier & Groth 2001, 2002; van Waerbeke et al. 2001a). 

Our method for addressing this problem makes use of the genus curves for weak lensing 
maps of the convergence field. Sato et al. (2001) have studied quantitatively the method 
based on two dimensional Minkowski functionals (the genus curve is one of the Minkowski 
functionals) via IV-body and ray-tracing simulations, and demonstrated its feasibility in the 
determination of the cosmological density parameter, Q m . Using the two dimensional genus 
statistic, Matsubara & Jain (2001) also studied the topology of weak lensing fields, and 
found that the genus curve for large smoothing angles agrees with the predictions from linear 
theory, while that for smoothing angles smaller than 10' shows the non-Gaussian signatures 
of gravitational clustering and differs for open and flat cold dark matter models. 

In the present paper, we report the first measurement of genus curves for the two- 
dimensional mass map in the neighborhood of distant rich galaxy cluster MS1054-03, recon¬ 
structed from weak lensing data obtained by Suprime-Cam (Subaru Prime Focus Camera) 
on 8.2m Subaru telescope (Urnetsu 2001; Umetsu et al. 2002). We examine the non- 
Gaussian feature in the measured genus curve by comparing it with a Gaussian prediction, 
and demonstrate the feasibility of the method to examine the statistical properties of the 
large scale mass distributions in the universe. 


2. Observation 

We observed the rich, X-ray-luminous galaxy cluster MS1054-03 at z = 0.83 on the 
night of 2001 March 29 using Suprime-Cam mounted at the prime focus of 8.2m Subaru 
telescope (Umetsu 2001; Umetsu et al. 2002). We obtained a mosaic of i? c -band images 
with a total exposure time of 4.6 ksec, which covers the sky area of ~ 30 x 30 aremin 2 
with the seeing FWHM of (P/8. For weak lensing analysis, we selected galaxies of R c < 26 
with signal-to-noise ratio (S/N) greater than 6 in the central 20 x 20 aremin 2 region. The 
number of objects in the weak-lensing catalog is 11609 (29.2 aremin 2 ). We have reconstructed 
the projected mass distribution from the image distortions of sample galaxies via Kaiser & 
Squires method (Kaiser & Squires 1993; Umetsu, Tada & Futamase 1999) in Fourier space. 
We first calculated the Gaussian smoothed shear field with d G = 0139 from the weak lensing 
catalog, where the Gaussian window is defined by Wq{6) = exp(— \9\ 2 /29q)/2tt9q. The shear 
map is then 100% zero-padded on each side to avoid spurious effects from periodic boundary 
condition. In Figure 1, we show the reconstructed convergence map with a side length of 19' 
and pixel size of 01074. The reconstructed convergence map shows a significant mass peak 
with S/N of 5.4 associated with the rich cluster MS1054-03. 
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3. Genus Statistic 

The genus statistic provides a topological description of large scale structures in the 
universe. Since the weak lensing field is mapped on the sky, the two dimensional genus 
statistic is used in this case. For the local convergence field k shown in Figure 1, we consider 
contour lines of a threshold value labelled by u. The threshold v is usually taken as n = iscr 0 
for random Gaussian field, where ay is the rms of the field, Uq = (k 2 ). The genus curve is 
expressed as function of this threshold level v. The two dimensional genus G 2 is intuitively 
defined as (Adler 1981; Coles 1988; Melott et al. 1989; Gott et al. 1990), 

G 2 = (number of contours surrounding regions 

higher than the threshold value) 

— (number of contours surrounding regions 

lower than the threshold value). (1) 


For a random Gaussian field, the genus curve is analytically given by (Matusbara 2000), 
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where ay is defined by erf = — (kV 2 k). 

Let us consider the meaning of value of the genus. If the value of genus is positive, then 
there is more isolated clusters than isolated voids and thus we can get an image that the 
topology is collection of isolated islands. So, when the value of genus takes large positive 
value, the topology is compared to ’meat-ball’. On the other hand, if the value of genus is 
negative, there is less isolated clusters than isolated voids. In this case the topology gives 
an image of the complicated shape with a lot of holes, and the topology is compared to 
’Swiss cheese’. Thus we say the genus curve shows a ’meat-ball shift’ and ’Swiss cheese shift’ 
when the curve deviates from Gaussian shape to positive direction and to negative direction, 
respectively. 

Gaussian expression of equation 2 is not expected for the present weak lensing data 
because the survey area is not large enough to probe the linear scales. However, it is known 
that the shape of the genus curve of mass distribution fields is not strongly affected by 
nonlinear gravitational evolution as shown by Gott, Weinberg & Melott (1987), Weinberg, 
Gott & Melott (1987) and Melott, Weinberg & Gott (1988). Thus we use the method of 
rescaled labeling of v. For non-Gaussian fields, the threshold is defined through the fraction 
of area / on the high-density side of the contour lines, which is related to the corresponding 



5 


value of u for a Gaussian field as, 

i r°° 

f — _ / e _<2 ^ 2 dt. (3) 

For a Gaussian field, this definition gives the same result as specifying the threshold directly 
in units of a 0 . On the other hand, for a non-Gaussian field, it compensates for the horizontal 
shift of the genus curve clue to the shift of distribution function of the density field. We 
adopt this prescription for setting the threshold values of k in the results shown below. 

To calculate the genus curves for the convergence map given as pixel data, we employed 
the method developed by Winitzki & Kosowsky (1997), Novikov, Feldman & Shandarin 
(1999) and Schmalzing et al. (1999). We confirmed that this code gives the same result 
as from CONTOUR2D developed by Mclott (1989). At each threshold level u, the genus 
is obtained as an average over three genus values with us shifted by 0 and ±0.1 as done in 
Park, Gott & Choi (2001). 


4. Results 

Figure 2 shows the genus curve measured from the convergence field shown in Figure 1. 
The Gaussian prediction curve is also calculated by equation 2, with the normalization factor 
obtained from the lcast-y 2 fitting to the measured genus curve. The errors for the genus 
measurements are estimated from 100-bootstrap resamplings of the background galaxies. 
We evaluated the curvature only at interior vertices for genus calculation. Figure 2 shows a 
clear difference between the random Gaussian prediction curve and the curve measured from 
the observational data. Furthermore, we found that compared with the random Gaussian 
prediction, the observed curve was shifted upwards near the threshold of u — 0. This 
indicates that the number of isolated clusters dominates over number of isolated voids, which 
is the so called ’meat-ball shift’. This phenomenon is usually seen when the field evolves 
into non-Gaussian from Gaussian. However, it is difficult to judge from this result shown in 
Figure 2 whether this non-Gaussianity for the genus curve is influenced by the rich cluster 
contained in this region or by the large scale structures. 

To clarify the origin of the non-Gaussian nature, we divided the reconstructed conver¬ 
gence map into four square regions with the same size, and calculated the genus curve for 
each sub-region. As indicated in Figure 3, we call each sub-region as region 1, 2, 3, and 4, 
respectively. The rich cluster MS1054-03 is contained in region 3. 

Figure 3 shows the genus curves for the sub-regions and their best-fit Gaussian genus 
curves, obtained in the same way as for the whole field. From Figure 3, we see that the genus 
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curve for the region 3 has a large skew and a ’meat-ball shift’. On the other hand, the genus 
curves for the regions 2 and 4 are consistent with random Gaussian curves. The genus curve 
for the region 1 shows a slight excess in the high threshold region. Our weak lensing analysis 
shows a significant mass concentration associated with the rich cluster, therefore the highly 
non-Gaussian feature in the genus curve could be induced by this rich cluster in the region 
3. 


The skewness is also a powerful statistical tool to measure the deviation from a random 
Gaussian distribution, defined by = (ft 3 ) / (ft 2 ) 2 . We summarized in Table 1 the skewness 
measured for each sub-region. We see the skewness for the region 3 is quite larger than that 
for any other region. 


5. Discussion and Conclusions 

We have measured for the first time the genus curves of the convergence map recon¬ 
structed from image distortions of faint background galaxies, as the first step for addressing 
the statistical properties of the mass distributions on large scales. The genus curve measured 
in the whole survey field (Figure 1) shows a deviation from the random Gaussian predic¬ 
tion, i.e. non-Gaussianity of the field. By dividing the whole map into four independent 
regions and analyzing them separately, we showed that the non-Gaussianity is induced by 
the rich cluster contained in the survey region (region 3). On the other hand, as for the 
regions without the rich cluster, we also found a slight excess of the genus amplitude in the 
high-thresholds for region 1, whereas the genus curves for regions 2 and 4 are consistent with 
random Gaussian distributions. We thus averaged the genus curves of the regions 1, 2, and 
4 by inverse variance weighting at each threshold bin. The resulting genus curve shown in 
Figure 4 is consistent with the random Gaussian one at each threshold bin within la error. 
We thus conclude the excess of the genus curve in region 1 is not due to the non-Gaussianity 
but is caused by the sample variance due to the small area of the region as discussed by Park 
et al. (2001). 

Our results might be influenced by the choice of smoothing angle ( 9q = 0f39 in the pa¬ 
per). The effective number of galaxies used for the local shear average decreases as smoothing 
angle decreases, and then the reconstructed convergence map is much influenced by the ran¬ 
dom noise due to the intrinsic distribution of galaxy cllipticities. Thus the genus curves 
approach to the random Gaussian predictions as smoothing angles become smaller. Matsub- 
ara & Jain (2001) discussed the effects of the smoothing angle on the genus measurements 
by the simulations with survey area of about 10 degree 2 . They concluded that for the 
smoothing angle of 9q = 1', the genus curve shows a non-Gaussian ’meat-ball shift’. On 
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smoothing angles smaller than 1', it is not clear yet whether the observed non-Gaussianity is 
real or not. We have not used 6q = V because it could decrease the number of independent 
samplings, and thus affect our analysis. 

In order to clearly see the deviation clue to large scale structure, we definitely need to 
observe much wider region which will be planned in the near future survey by Suprime-Cam 
on Subaru telescope. The present study demonstrates the feasibility of the genus statistic 
for the weak lensing mass maps to examine directly the statistical properties of the large 
scale mass distributions in the universe. 
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pixel value 


Fig. 1.— The weak lensing convergence map reconstructed from the image distortions of 
background galaxies in the i? c -band image of Suprime-Cam with Gaussian smoothing scale 
of 0f39. The side length is 19' and pixel size is 0^074. The four sub-regions are also indicated 
in the figure. The rich galaxy cluster is contained in the region 3, indicated by cross. 
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Fig. 2.— Genus curve per steradian of the convergence field shown in Figure 1. The open 
squares and dashed-curve indicate the genus measurements from the observational data. 
The solid curve indicates the random Gaussian prediction with best-fit normalization to the 
measured genus curve. 




Same as Figure 2 but for each of the four sub-regions indicated in Figure 1 




o 



V 

Fig. 4.— Same as Figure 2 but for the region without the rich galaxy cluster MS1054- 
03, obtained by averaging the genus curves for the regions 1, 2, and 4 by inverse variance 
weighting at each threshold bin. The genus curve is consistent with the random Gaussian 
prediction. 
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whole field and the four sub- 
100 bootstrap resamplings of 

galaxies. 


region 

skewness 

whole 

8.53 ± 5.99 

1 

4.68 ± 7.42 

2 

1.62 ± 11.32 

3 

19.45 ± 10.02 

4 

2.42 ± 10.76 


Table 1: Skewness of the convergence field measured for the 
regions shown in Figure 1. The errors are estimated from 



